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CTS \ We develop a coordinate space renormalization of massless Quantum Electrodynamics using the 

T^ ' powerful method of differential renormalization. Bare one-loop amplitudes are finite at non-coincident 

external points, but do not accept a Fourier transform into momentum space. The method provides 
a systematic procedure to obtain one-loop renormalized amplitudes with finite Fourier transforms in 

,^ ■ strictly four dimensions without the appearance of integrals or the use of a regulator. Higher loops 

are solved similarly by renormalizing from the inner singularities outwards to the global one. We 

k> , compute all 1- and 2-loop IP! diagrams, run renormalization group equations on them and check 

H \ Ward identities. The method furthermore allows us to discern a particular pattern of renormalization 

under which certain amplitudes are seen not to contain higher-loop leading logarithms. We finally 
present the computation of the chiral triangle showing that differential renormalization emerges as a 
natural scheme to tackle 75 problems. 
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0. Introduction 

All one-loop diagrams in Quantum Electrodynamics (QED) are products of propagators when 
written in coordinate space and, therefore, finite for separate external points. Apart from tadpoles, 
this is in fact true for any quantum field theory. This striking observation has been in the mind of 
the physics community since long ago though little practical use of it has been made. The need for 
renormalization comes through the fact that the product of propagators is not a distribution and, 
thus, has no Fourier transform. In crude terms, renormalization is a procedure to extend products of 
distributions into distributions. 

Differential Renormalization (DR) [1] provides a prescription to implement this project which 
is in a sense minimal. The basic idea is to write an amplitude in position space and then express 
it as derivatives of a distribution less divergent at coincidence points. The typical computation of 
primitively divergent Feynman integrals is substituted for solving trivial differential equations. This 
prescription is complemented with the instruction of using the derivative in front of the distribution 
as acting to its left, which is standard in the theory of distributions. 

We claim that differential renormalization is minimal because it never changes the value of a 
primitively divergent Feynman diagram away from its singularities, and neither does it modify the 
dimensionality of space-time or introduce a regulator. It is simple and addresses straighforwardly the 
concept of renormalization, rather than regularization. In spirit, it is close to BPHZ renormalization 
[2], in that both prescriptions deliver renormalized amplitudes graph by graph. In the absence of a 
renormalized lagrangian, DR needs a supplementary effort to prove unitarity. In [3], this problem has 
been analyzed for A(/)''-theory and perturbative unitarity proven to 3-loop order. 

A number of theories have been used as a test for differential renormalization. In [4], the 
supersymmetric Wess-Zumino model was studied up to three loops with ease. Conformal invariance 
in QCD was exploited in [5] thanks to the position space nature of differential renormalization. The 
method also extends to massive theories [6] and to lower dimensional theories [7]. 

Other coordinate space regularization and renormalization methods have occasionally been used in 
the past. From Schwinger [8] to all the recent conformal field theory developments [9], the advantages 
of postponing loop integration to higher loops and the possibility of exploiting conformal invariance 
have led to numerous applications of coordinate space methods. 



In this paper we present a comprehensive study of the differential renormalization of massless 
QED up to two loops. There are obvious reasons to undertake such a computation. QED is a physical 
theory with gauge symmetry. The consistency of the method passes a stringent test because of the 
existence of Ward identities (Wis) and potential problems hidden in overlapping divergences. It also 
faces the problem of the chiral anomaly and the correct treatment of 75. 

Our results are in perfect agreement with other computational schemes (e.g., [8], [10], [11], [12], [13]) 
but are obtained in a simpler form. We devote Sections 1 and 2 to presenting the renormalized am- 
plitudes up to two loops. Whereas several standard treatments of QED will perform the subtractions 
on mass-shell, here renormalized amplitudes are subtracted at an arbitrary and unphysical scale M, 
due to the fact that we are dealing with massless QED; this massless limit exists as a quantum field 
theory and provides substantial information about massive QED. We pay special attention to the way 
differential renormalization uses the basic Wis of the theory and how they are repeatedly checked in 
two-loop computations. Deeply rooted in the nature of DR lies the idea that each Feynman diagram 
has as many hidden scales as singular regions. Some of these scales are fixed by Wis. The rest are 
renormalization-scheme prescriptions. Keeping all this freedom manifest allows the method to "pre- 
dict" that the renormalization group equations will display a particular organization. This "structured 
renormalization group", discussed in Section 3, provides a deep insight into the absence of promotion 
of leading logs in some diagrams, which is seen as accidental in standard treatments. Finally, we 
review the calculation of the chiral anomaly in our regulator-free way in Section 5. We also compute 
the absence of infinite renormalizations (i.e., absence of In's) at two loops of the triangle amplitude. 
As differential renormalization never leaves four dimensions, it stands as a good candidate to work 
with 75-related observables. The anomaly is just one instance. 

1. One-loop Renormalization 

1.1 Conventions and definitions 

\n d ~ 4 Euclidean space, massless QED is described by the Lagrangian: 

^ ^ i^''^ + ^(a^A^)2 +^(^ + *e4)V', (1-1) 



where j4^(a;) and tl^{x) are the usual C/(l) gauge and Fermi fields, F^,^{x) = dnAiy{x) — di,Afj_{x) 
is the gauge field strength, a is a gauge fixing parameter, and 7-matrices satisfy the usual Clifford 
algebra, {7^,7!/} = 2(5^i.. 

Loop diagrams are calculated with the following Feynman rules: 
i) the gauge field and massless fermion propagators are, respectively. 



G|.°.n--y;«) = 4^[( 



l + a 



2 J (x-y) 



V _^^^_dx-y)t^ix-y)-^ 



(x - yy 



^^°^(--y) = -i^7.5,^^^, 



ii) to each vertex is associated a factor ie7^, 

in) internal coordinates are integrated over, 

iv) closed fermion loops are multiplied by a factor —1, and 

v) diagrams have no symmetry factors. 

We also use the convention that for derivative operators O = 9^,^,n, etc., we take Of{x — y) to 
mean 0(x)f{x — y), unless explicitly stated otherwise. 

Renormalization of QED entails the renormalization of only three different vertex functions: the 
2-point vacuum polarization, n^^(x — y), the 2-point fermion self-energy, Y,{x — y), and the 3-point 
vertex, V^{x — z,y — z). These will contribute to the effective action in the following way: 

T,s = j d\d\j\^A^{x) (Al-^)d^d,-5^,n\5^^\x-y)-n^,{x-y) A,{y) 

+ ^(x) (^,5W(x -y)- ^{x - y)) ^(y)} 
+ f d^xd^d^i^ix) hej^S(^\x - y)S<^^'>ix - z) + V^{x - z,y - z)] ^{y)A^,{z) + ■■■. (1.2) 



As defined above, these renormalization parts will contain only IP! loop contributions and II^j^ 
and S will, furthermore, lead to the following full photon and electron propagators: 

s = ^(0) + 5(0) . s . 5(") + 5(") ■ s • ^w • s • 5(") + • • • = (5(")-i - s)-\ 



(1.3) 



where the dot indicates a convolution. To all loop orders, renormalized vertex parts satisfy renor- 
malization group (RG) equations [14] and, as a consequence of gauge invariance, the following Ward 
identities: 

d^ll^,{x-y)^0, (1.4) 



dzf^ 



■£^V^{x -z,y-z) = -te[6^^\z - x) - 5^^^\z - y)]S(x - y). (1.5) 



The specific coefficients in the second Ward identity are uniquely fixed by the tree-level identity solved 

by yj"' = «e7^(5(4)(a; - y)5^^){x - z) and S^^^-^ = ^ S^^^x - y). 

In what follows, we shall find the 1- and 2-loop corrections to H^^, S and V^. We will verify the 
above Ward identities, and the requirement that renormalization parts satisfy renormalization group 
equations will allow us to calculate the renormalization group functions of QED (beta-functions and 
anomalous dimensions). 

1.2 Renormalization 

The bare 1-loop vacuum polarization (Fig. la) reads: 

' ' 1 .1 



^T'^^'^ix -y) = -(^] Tr (7,^ 7-—T,lu9 



An^ J \ [x — yY (2/^2;) 



(1.6) 



with a = e^/An the fine structure constant. Using the basic DR identity (cf. Appendix A for all the 
identities needed in renormalizing amplitudes up to two loops): 



1 _ 1 _ Inx^M^ 



- = --n —^ , (1.7) 



we find the renormalized value of the one-loop vacuum polarization: 

^(- - V) - ^(d.d. - s,.nW-^^^^^^ , (1.8) 

with the proviso that the total derivatives above, as well as in all other renormalized amplitudes 
throughout, are to be understood as acting to the left. Here and below we adopt the convention 
of appending a subscript (H, E or V) to the renormalization scales appearing in the different renor- 
malization parts, since these scales are a priori independent (but eventually related through Ward 
identities, as we shall see). This quantity gives us then the renormalized AA 2-point function: 

rt^i^ -y) = ((1 - \)d,d. Vn) 5^'\^ -v)~ nW(x - y). (1.9) 



We now impose that it satisfy the usual RG equation: 

where /3(a) is the QED /3-function, 7^(0) is the anomalous dimension of the gauge field, and 7a(a) 
a "/3-function" associated to the running of the gauge parameter a. This is then easily seen to lead 

to the 1-loop values: 

a 

^A{a) = — 

and 

7a(«)--^ (1-11) 

{I3{a) has dropped out above because it is of higher order in a). 

It is well-known that the above computation can be used in the framework of the background 
field method to find the beta-function of QED [15]. It turns out that the relation P{a) = 2aj{a) 
holds to all orders. Therefore, we can anticipate that the beta function at one-loop is 

/3(a) = ^. (1.12) 

This result will later be confirmed by an independent computation of the one-loop vertex. 
The bare 1-loop fermion self-energy (Fig. lb) is: 



,u..-.,. -„ = (_£,)%„ ^i,„ 



1 + a Sfj^iy ^-^_,{x-y)f,(x-y)^ 



(1.13) 



2 (a; — y)^ {x — y)'^ 

The renormalization of this amplitude is straightforward, apart from the following important subtlety: 
if we use the same mass scale A/s in renormalizing the two parts coming from the two different pieces 
in the photon propagator, we will end up with a renormalized amplitude which is directly proportional 
to the gauge parameter a; now, in Landau gauge (a = 0), it would then vanish entirely, and it is easy 
to see this would be inconsistent with the Ward identity. It is in fact a known property of QED that, 
in Landau gauge, although the bare 1-loop fermion self-energy vanishes, the renormalized self-energy 
does not. In momentum space, it is equal to a finite constant times ^, and this is a reflection of 
the ambiguity inherent in the linearly divergent integral defining the amplitude. This turns out to be 
transparent in differential renormalization: in renormalizing the two pieces in question, no one tells 
us we should take the same mass scale /Vfs, and so we do not. We renormalize the first piece (the 



Feynman gauge part) with a scale Ms, and the second piece with a scale M^ related to the first one 
by InAfg = InM^ — A. This then leads to the following renormalized 1-loop fermion self-energy: 



IGtt'^ \ (x — y) 



2 ..2 . . . (1-14) 



«. ln(x-,)M|_Aa(l-a)^^(,)^^_^^ 



IGtt'^ (x — y)2 Stt 

As expected, in Landau gauge, all that survives is precisely the contact term mentioned above. We 
will see shortly that the Ward identity will not only determine a relation between the self-energy and 
vertex mass scales, but also a unique value for A. The renormalized ^''0 2-point function, 

V^'^ix - y) = ^3^'^\x - y) - T.^^Hx - y), (1.15) 

then satisfies the RG equation: 

(^^a^ + ^^"^^ + ^"^"^^ " ^^'^^"O ^'^^^^ - y) = 0' (1-16) 

and this will lead to the following 1-loop fermion anomalous dimension: 

, . act /-, -,-N 

7^(a) = -^ (1-17) 

(here, both f3 and ja drop out). 

Finally, we now renormalize the 1-loop vertex (Fig. Ic). Its bare expression is: 

V;«'^-°(x-z,y-z) = 



I (7p7a7p7b7<T)9a7 77 4- 



47rV '"""'^"""' "(x-z)2^''(z-y) 



1 + a^ Spcr _^ {x~y)p{x-y). 



2 / (a; - y)^ {x - y) 
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2 /Cl\3/2 11 no 

= -^^ \-) ilplalplbla) 9a , _^p 9b . _ ^ [(« - l)9p9ff + (Jp^o] ln(x - y) fi . 

(1.18) 

Like for the fermion self-energy, the photon propagator in a generic gauge leads to two pieces to 
be renormalized, in principle with mass scales My and My. However, it turns out that this need not be 
done here, that is, the choice of the same My throughout will not lead to any inconsistencies (still, one 
may choose different scales at will; this would simply correspond to different choices of scheme). We 
briefly sketch the procedure involved in renormalizing the expression above. In Feynman gauge (a = 1) 
it is most transparent [1]: there are three propagators l/{x — x')'^ forming a triangle, with derivatives 
acting on two of them. This has dimension L^^ and therefore it is log divergent in the ultraviolet. 
One first integrates by parts until the two derivatives are acting on the same leg, say, dadt ,^} y^ , 
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and then subtracts and adds a trace piece thus: (dadb — jSabn) , ^ ^ + \5ab^ i 1 -,2 - The surface 

4 [^x y ) 4 ^2; y) 

terms, with total derivatives acting outside, is power counting L^'' and thus finite, and the traceless 
combination of derivatives is also finite (due to tracelessness). The divergence has been isolated in 



the term j5abU-r-^. 



(x-y)^ 



7r^(5ab(5^^'(x — y), which turns out to be easily renormalizable through 



the DR identity, Eq.(1.7). For a generic gauge, the same principle of separating divergent terms into 
trace and traceless pieces applies, only in this case the 7-matrix structure in front complicates things 
a little: we integrate da and db above by parts onto the photon leg, and then look for the coefficient 
A to make the expression 



ilplaliilbla) {dadb + MabU) [{a - Vjdpda + 5 p„u\ 



(1.19) 



a traceless (and thus finite) combination of derivatives. The appropriate value is ^ = —a/{a + 3); 
adding and subtracting that trace piece, we are able to find the renormalized expression for the 1-loop 
vertex: 

V<^'\x-z,y-z) = 



dy 
16 



1 



db- 



1 



[x-zY {z-y) 



[{a - l)dpda + Span] ln(a; - yf^i 



+ 4 [(a - 1)7^ 7h - 27h7^] ^_ ^ 



1 



1 



(x-z)2(z-y): 



r7<T 



1 



dpdcr - jSf_,an 



{x — zY{z — yY {x — yY 

1 . 2 \n{x-zfM^ 



[x-yf 



^TT a-ffU - 



(x — zY 



5^'\x-y) 



(1.20) 



We remark here on the fact that the renormalization piece above (containing InA/y) is directly 
proportional to the gauge parameter a, and this is also true of YJ^^^{x — y) found above. This is a 
reflection of the well-known fact that, apart from vacuum polarization infinities, QED is 1-loop finite 
in Landau gauge. 

The renormalized 3-point vertex 

Vf^'^ix ~z,y-z) = ie^p6^^Hx - y)S^^Hx - z) + V^^'^x -z,y-z) (1.21) 

satisfies the RG equation: 



^^My + ^^"^^ + ^'^^"^^ " ^^'^^"^ " ^-^^"0 ^'^'''^^^ ~^,y-z)^o 



(1.22) 



with values found previously for 7^ and ja, and 



/3(a) 



2a' 
Sir 



(1.23) 



(and, again, 7a is of higher order than we are considering). This confirms the result from the 
background field method shown previously. Naturally, in the above equation the a-dependent pieces 
(Mg^ and 7^) and the a-independent pieces {li-§^ and 7^1) cancel separately. 
We now consider the Ward identity for V^ and S. 

1.3 Ward Identity 

For separate points x ^ y ^ z, Eq.(1.5) is trivially verified. The subtlety involved in a coordinate 
space Ward identity is, however, its validity at contact. In order to verify Eq.(1.5), then, we must be 
careful in particular not to lose any contact terms due to formal manipulations with delta functions, 
and the best way to guarantee a correct procedure is to integrate it over one of the external variables. 
This was done in [1] for Feynman gauge, and we use the same procedure here for an arbitrary gauge. 
The integrated form of the Ward identity is; 

j d\j^V^^\x -z,y-z)= le^ix - z). (1.24) 

Integrating the expression for the renormalized vertex, Eq.(1.20), over y, one finds: 

l + 2a \ 1 a h\{x- zfM'^'^ 



Av«(.-...-=)^j;i(f)"^[(<J,^-l±^.„n) 



[x-zf ' 2"^"°" {x-zY 

(1.25) 



Acting now with the z derivative: 



^^^nl^n--,y-) = ^^D(^ — \:;r^ — '-) ■ (1-26) 

Now, we compare this with: 

.eS(^)(. -z)-^jn ^ ^H:c-zfMl^Kl-a)/2 ^ (,.27) 

Setting a = gives A = 3, and then, for any a, we also find: 

Ml I , , 

lnT7T = ;^- 1-28 

M2 2 ^ ' 

We see then that a Ward identity relating two renormalized amplitudes will, in the context of 

differential renormalization, enforce a relation between the scales that renormalize the amplitudes. 

It is an important point we will analyze further that the converse is also true, namely, mass scales 

that renormalize amplitudes not related by a symmetry are not related, and this in turn also enforces 

constraints in the renormalization of these amplitudes. At two loops, in particular for the anomalous 

triangle, the use of the above mass relation will be crucial for the consistency of our calculations. 
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2. Two-loop Renormalization 

At two loops, the IPI diagrams renormalizing 11^,^, E and V^ are depicted in Fig. 2. Throughout 
this section, we will perform all calculations in Feynman gauge (a = 1); since 7a^ acting on any 
2-loop diagram will be one order higher in a, this will not affect the verification of RG equations. 

2.1 Vacuum Polarization 

We begin with the simplest of the two vacuum polarization diagrams (Fig. 2a), and from now 
on we will use translation invariance to set one external point to zero in all diagrams, for the sake of 
simplicity. Its bare expression is: 

<-)--(.) -0^Jd\d% Tr (7.^ (^S^^^(«--)^^>^i) . (2.1) 

Standard manipulations lead to: 



n(f;)'^-<=(x) 



e 



4 



487r6 



Vd.0.-s,.n)^'^^'''--'/'^ '' 



^V-P-. .^.-y, ^4 , ^^ 



(2.2) 



We note that this diagram is not transverse by itself; the non-transverse piece will be cancelled 
when we consider the entire 2-loop vacuum polarization. The renormalization now proceeds with the 
straightforward use of the DR identities listed in Appendix A, and we find: 






(2.3) 
where M is a new renormalization mass parameter appearing at two loops. The diagram with the 
fermion self-energy inserted on the lower leg of the loop will have the same value as this one. We 
note here that the lnM| coming from the 1-loop self-energy subdivergence has been promoted to a 
In Mg at two loops: this is the self-consistency of the renormalization group at work. We will choose 
everywhere the same 2-loop mass scale M. In renormalization terms, this simply corresponds to some 
choice of renormalization scheme: because these appear as D '" ^ /^ , different values of M will lead 
to amplitudes that differ by finite contact terms or, in other words, by finite renormalizations (and 
this will be true for the new M's appearing at each loop order). Of course, once we set the (2-loop) 
M's we like for the renormalization of S and V^ in particular, they will subsequently be related by 
the Ward identity (but this will not be of concern to us here, as we will not compute the 2-loop mass 
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relation). Throughout this section, then, we will only care to distinguish between Ms, My or Mu 
and other mass scales M when the former appear as promoted In's, i.e., as In^ A'l^, etc.. 

The other vacuum polarization diagram. Fig. 2b, is the most difficult integral we have had to 
perform. The bare expression is: 



H- i^/^jj.z^i> J y f (x — u)^ [u — y)^ [v — uY [v — xY J yu — v)'^ 

(2.4) 

There are potential problems related to overlapping divergences in this amplitude, and so we must 
examine how differential renormalization deals with them. There are two subdiagrams which contain 
log singularities related to the regions u ^ v ^ x and u ^ v ^ Q. We remove both divergences by 
pulling out derivatives in x and y in a symmetric way as we did in the case of the one-loop vertex, 
Eq.(1.20). More explicitly, the important intermediate step is 

"d-, T^Oa-, ^Ob-, ZT^Oc 



2 



' {x — uY [y — uY {y — vY {x — vY {u — v) 

. d 1 1 ^ (f, P, ^"^"^n) ^ 2e S'^JX-V) , 1 

^dxd\x-uY "(x-vY' {x-uY^"" 4°'{x-vY "Ux-uY^iu-vY ^ ' 

r ^ ^ 1 a ^ ^ ^ ^a a _ ^n^^^^ ^ 2r ^Hy_^. 

^dyb^y-vY {y-uY' {y-vY^ 4°' {y - uY {y - vY ^ ' 

This manipulation makes it evident that each subdivergence is cured separately thanks to the fact that 

in coordinate space the external points are kept apart (whereas, in momentum space, the momentum 

integral would make the two singularities overlap). The rest of the computation, carried out with the 

help of Feynman parameters, although somewhat lengthy, is conceptually simple since only a global 

divergence needs further correction. The final renormalized value is: 



(2.6) 

Taking into account the mass relation (Ward identity), Eq.(1.28), the entire 2-loop renormalized 
vacuum polarization reads: 

n(i)(x) = 2n%^){x) + n(i^)(.) = ^(^)'(a,9. - s,^n)n^^^ . (2.7) 

This is automatically transverse, and furthermore the In contributions have cancelled. This is an 
important point which we will further elaborate in Section 3. 
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The 2-loop renormalized AA 2-point function 

r;!.^W = f (1 - hd,d. 6,^) ,5W(x) - n«(x) - u^^}{x) (2.8) 



will satisfy an RG equation, Eq.(l.lO), with Mn g^ substituted forMn^^^ — ^-^WT ('" g^n^i'sl. 3*^ 
RG equation will include the sum of the derivatives w.r.t. all masses present in an amplitude). This 
will yield the 2-loop RG functions: 

2 

a a 



and 



^^(") = 3^ + 4^2 



7a(a) = -^-^ = -2a7A(a). (2.9) 



As mentioned before, the same calculation in the background field method would have led us to the 
2-loop beta function, which we will derive and present independently through the calculation of the 
2-loop vertex in Sec. 2.3. 

2.2 Fermion Self-Energy 

Next, we consider 2-loop contributions to the fermion self-energy, shown in Figs. 2c-2e. In order 
to eventually verify RG equations and Ward identities, we also give the logarithmic mass derivative 
Mg|j acting on each one of the amplitudes. The first diagram is that of Fig. 2c: 

S*"^'^^°(-) = -(^/^'-'^'-^7.^(^^sW(--^')^^7m. (2.10) 

Again, integration by parts, properties of 7-matrices and DR identities lead to the renormalized 
value: 

^""'-'-T^e'^°( '°'^'"t:'"'"> 

where we have also used the identity 

'a. ftf!M) . I9. f '"'^1-'/^ ) . ,2.12) 



The mass derivative of this amplitude is: 

dM 27r 167r2 ^ ' 

The next contribution is that of Fig. 2d: 

j.(2.)bare(^) = - y^ / d^ A 1 ^^7p^ ^^7.^ ^^7^^ 4> ' (2-14) 

(47r^)-^ J u^ (x — vY [x — uY (u — vY v 
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This amplitude is solved by using vertex-like manipulations. It turns out that no new integrals are 
needed apart from those appearing in the two-loop vacuum polarization diagram 2b. We also make 
use of the following trick 



-An' ^ '- + H^ + 2da- -^da 



which simplifies part of the computation. The renormalized value we finally get is: 

>2„ hi^x'Ml 



and the mass derivative of this amplitude is: 

M-^s(2'i) = -E(i)-^S(0), (2.16) 

dM TT Sir' ^ ^ 

where we have used the 1-loop mass relation to express S'^^ with the mass scale Afs rather than 
Mv 

We finally present the last of the 2-loop fermion self-energy diagrams (Fig. 2e): 

^*"^'^^^^(-) = -(S5yj7.^^>/rf^rf^(^n(VH.-«)l . (2.17) 

The renormalization of this amplitude is straightforward, and we only point out that while there is in 

(1) 2 

principle the possibility that the InAfn in 11)^^: would get promoted to a In Mn, the 7-matrices in 

(1) 2 

the amplitude, together with the transverse operator in nj^J, conspire to simply cancel all In 's. This 

(2) 
is again an instance of the feature we have seen previously in XI/i^;, namely, the apparently unexpected 

cancellation of some particular divergences. The final, renormalized result is: 

^^^^H^)--^,{-)^u'^^. (2.18) 

■Utt'' \7r/ X'' 

The mass derivative of this amplitude will be: 

dM 47r2 ^ ' 

Whereas we would generally expect the mass derivative of a 2-loop amplitude to generate 1-loop and 
tree-level amplitudes, we see this does not happen here. In Section 3, we will examine this more 
closely. 
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We now add all these amplitudes to find the 2-loop renormalized ifjip 2-point function: 

a fa\ ^ Inx'^Ml 1 /a\2^ fln^x'^Ml-Tlnx'^AP 

)--^r^[-)^n — ^-TT^^T - ^° ^ 

Idtt"^ Vtt/ x^ 1287r^ Vtt/ \ x^ 

(2.20) 

We have the 1-loop correction in a generic gauge a, and the 2-loop terms in Feynman gauge. This 

is written thus because in RG equations we will need ^ on the 1-loop term but not on the 2-loop 

terms. The 2-loop fermion anomalous dimension given by the RG equations is (in Feynman gauge); 



a 3q:^ 
4^ ^ 32^2 



7^(") = — -:;TT- (2-21) 



2.3 Vertex 

We now turn to the computation of 2-loop corrections to the gauge coupling vertex, shown in 
Figs. 2f-l. Many of the integrals are extremely difficult to perform analytically, but fortunately they 
need not be done for the purposes of verifying RG equations and Ward identities (although we will 
not present the 2-loop relation between My and M^). The parts of these diagrams that do need 
to be computed fully are the renormalization pieces from both global and internal divergences, and 
those are reasonably simple to calculate. Because the expressions for the renormalized amplitudes are 
somewhat lengthy, we present them in Appendix B. Here, instead, we will explain briefly the different 
techniques we have used and subtleties involved as we go along. Again, for the purpose of later 
verifying RG equations and Ward identities, we also present here the result of the logarithmic mass 
derivative A^ofy acting on each one of the renormalized amplitudes. We set z = in all amplitudes. 

The first diagram we consider is Fig. 2f: 

''^''''''''^^^y^--WYj^^ j '"^'"^^^^h^^^ (2.22) 

The standard integrations by parts and separation into trace and traceless pieces which were used 
to renormalize the 1-loop vertex are used here as well, without any added complication. We also need 
to consider the diagram identical to the one above, but with the self-energy insertion on the other 
fermion leg. Fig. 2g. The procedure is identical to the one used above, and its contribution to the 
RG equations at two loops is also the same. From the renormalized expression given in Appendix B, 
we can calculate the mass derivative of this amplitude: 
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Here and below, we are always taking Vp, in Feynman gauge, unless otherwise explicitly stated; we 
must also be careful to express the result in terms of My and not Afg (in the above, for instance, 
neglecting this would lead to a different coefficient for the contact piece Vp, ). 

The second vertex diagram we compute. Fig. 2h, contains a vacuum polarization insertion: 

^^(2.)ba.c(^^ ,) = ^ (7p7a7.7.7.) 9.^a,l j d\,d% (^^nW (u .)^-^ . (2.24) 

To renormalize this, we integrate the derivatives da and db by parts onto the photon leg. The transver- 
sality of njjj and the gamma structure in front will arrange things so that the resulting combination 
of derivatives will automatically be traceless, thus obviating the need for further renormalization and 
(again) avoiding the promotion of In Afn- This in turn leads to a peculiar form for the mass derivative; 

Now, we are left with the more difficult diagrams, which cannot be computed in closed form. 
Starting with the diagram of Fig. 2i, its bare value is; 

— I dSldSj- TiJ^aT -TT^db^rdc^rdd-, nj . 

(x — yy J v"- ^ ^) [^ ^ u) "^ "^ ['^^ ^ y) 

We first of all renormalize the upper vertex (connecting points z = 0,u, and f), by integrating by 
parts in z and separating the term integrated by parts into a trace and traceless piece (in the indices 
b and c). The surface term is finite by power counting, and the trace piece easily renormalizes to a 
structure of the form 

{ldltj,la) 7 yT da^ dd 2—^ , (2.27) 



which needs to be renormalized one more time (integration by parts and separation into trace and 
traceless pieces), giving a In^ promotion through the use of a DR identity. This is easily done, and 
now we are finally left with the global divergence {x ^ y ^ 0) present in the traceless piece (in be) 
we got after the first integration by parts. That term is power counting log divergent, and one might 
think the traceless combination of derivatives then makes it finite, but this is not so: the point is 
that there are other free indices in that expression (a and d). The integration is indeed made finite 
by tracelessness in the be indices, because the a and d derivatives can be brought out of the integral, 
but there remains a global divergence as a; ^ y ~ 0, because the expression is not traceless in all 
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indices. Attempting to separate the trace pieces from four free indices is hopeless because we cannot 
even do the integral in u and v, and so we resort to a technique used to renormalize the nonplanar 
3-loop 4-point diagram in Xcf)^, valid for primitively divergent expressions in general. Details can be 
found in [1], and we do not give them here. The idea consists in writing the factor ,^ -.^ in front as 

-{x-y) \^ ..^4 - . a V^ .^. • (2-28) 



{x — y)^ ' {x — y)^ 4 [^ ^ nY 

One can verify that the simple substitution above, as is, suffices to renormalize the global di- 
vergence we had (that is, that term will have a well-defined Fourier transform). The only subtlety 
in applying the mass derivative to the resulting renormalized amplitude lies in fact in this globally 
divergent piece. M-Mj on the term above is proportional to 5'^^^{x — y) and, like in [1], one must verify 
that what multiplies this is a representation of 5^^''{y — z) when x — > y. Adding up all contributions, 
the mass derivative then reads: 

M— T/(2') - — y(i) + -^ y(") (2 291 

dM > " 27r > + 16^2 V • l^-^^i 

The next diagram is that of Fig. 2j, with bare value: 

^5 



— fie) 

9a^T d\id\)db^T-. TT^dc- —dd- 

x'^ J W' [X ~ V) 



(2.30) 



(x — v)^ (u — v)^ {v — y)^ {u — yY 

The first step in renormalizing this diagram is integration by parts of the dd and dc derivatives, 
and separation into trace and traceless parts. There will be two of each; the first trace piece is easily 
renormalizable with standard DR identities, and the second one involves the structure 

{-fdlalii)da^dd^K{x,y), (2.31) 

where 

^^^^y^-j'^'^ ix-uYly-uYu^ - ^'-''^ 

This is renormalized one more time by integrating the derivatives by parts onto the X-function and 
separating trace and traceless pieces. For the first time, we encounter renormalization structures that 
do not correspond to any lower loop diagrams, and because of consistency with RG equations, these 
must cancel in the only other diagram left to compute. Fig. 21 (the nonplanar diagram). We find that 
indeed they do. The only other divergent pieces remaining are the ones traceless in cb and cd, which 
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can have a global divergence for the same reason as the equivalent term in the previous diagram: 
the presence of extra indices. These are primitively divergent, like in the previous diagram, and are 
solved in exactly the same way. Naturally, the diagram identical to this one, but with the vertex 
subdivergence on the opposite fermion leg. Fig. 2k, is solved in the same fashion, and leads to the 
same contribution to the 2-loop RG equations. We present the result of the mass derivative on these 
diagrams after the inclusion of the nonplanar diagram. Fig. 21, precisely because there are structures 
in these three diagrams which do not correspond to any lower loop diagram, and which cancel when 
they are added. So, finally, the nonplanar diagram is: 

(2.33) 






dSid V da - 






' [x — uY u^ v^ {v — y)^ (x — v)^ (u — y)^ 
To renormalize this, we integrate dc by parts around z = 0, and again separate trace and traceless 
pieces. Within the trace pieces, we will find the structure 



7m 



d 



dx'^ yx^y"^ dy°- 






x^y"^ dx'^ dy"^ 



(2.34) 



and these will precisely cancel with structures found in the two previous diagrams. We can now present 
the mass derivative acting on the sum of these three diagrams: 



d 



M-— [y(2j) + y(2fe) _^ y(20l ^ !1 y(l) _ 1±_ yiO) 



dM 



Stt 



(2.35) 



We are now ready to consider RG equations. These are easy to verify once we use the mass 
derivatives given above for the 2-loop vertices. The 2-loop Y'j^^^ vertex function, 

satisfies an RG equation, Eq.(1.22), confirming all the results given previously for different RG func- 
tions, and yielding also the 2-loop /3-function: 



/3(a) 



2a^ Qf^ 

"3^^ 2^' 



(2.37) 



This matches the value gotten by a background field calculation on the 2-loop vacuum polarization. 



2.4 Ward identity 
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At two loops, the method employed previously to verify Ward identities, viz., integrating over an 
external variable, becomes computationally difficult and not very illuminating. For our purposes here, 
we shall consider instead the following simpler procedure: we apply the logarithmic mass derivative 
to both sides of the (2-loop) Ward identity; as we have seen above, this yields 1-loop and tree-level 
vertices and self-energies, whose Ward identities have already been verified, so that our problem is 
reduced to that of matching the coefficients of 1-loop and tree-level quantities on both sides of the 
identity. We have, on the one hand, gained much in simplicity but, of course, on the other hand, 
some information contained in the Ward identity will thus be lost, namely, the contact terms coming 
from the finite parts of vertices, since these latter have no mass scales and will vanish when the mass 
derivative is applied. Although in this way we cannot derive mass relations like we did at one loop, 
the match we find represents a highly nontrivial consistency check of our 2-loop computations. 

Adding up all the contributions from Eqs.(2.13), (2.16) and (2.19), we find: 

^d ^(2c+2d+2e) ^ ^ ^(1) _ 2^ j.(0) (2.38) 

dM 27r 167r2 ^ ' 

and from Eqs.(2.23), (2.25), (2.29), (2.35), 

M— T/(2/+2s+2h+2j+2j+2fe+20 _ _^ y(l) _ ^Q ^(0) _ ^T/(l)/„ _ r)^ (o oq^ 



Once we use the fact that, from Eq.(1.25) 



(2.40) 



we immediately find the amplitudes we have calculated above do indeed satisfy the 2-loop Ward 
identity. It is also worthwhile noting, furthermore, that in fact we can divide the above 2-loop 
amplitudes into three sets that separately verify the Ward identity. They are: 

d 



d 
d 



V^'f+'^+'^\x,y) ^ -^e[S('Hx) S^^\y)]J:('^\x y) , (2.41) 

^(2,+2fe+2/) (^^ y^ ^ _^g[^(4) (^) _ ^(4) (y)]s(2rf) (x-y), (2.42) 



g^X''^{x,y) = -ze[d^^\x) - 5^^\y)]J:^'^\x - y) . (2.43) 

The vertex diagrams in each one of these sets are generated by attaching an external photon line 
in every possible way to an internal fermion line of the corresponding self-energy diagram. This is a 
vestige of the fact that each bare vertex - individually - formally satisfies a Ward identity with the 
self-energy gotten by eliminating the external photon line from that vertex diagram. 
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3. Structured Renormalization Group 

In this section we point out a certain structure exhibited by the renormalization of the different 
relevant vertex functions of QED. For that purpose, we gather here the results of the mass derivative 
of all 2-loop amplitudes. 

Vacuum polarization: 

mA [2n(,t") + n(i^)] = -$^{d,d,. 5,.u) Six) (3.1) 



Self-energy: 



dM ' "" f"" ' 27r2 



dM 2tt 167r2 "- ^ 

^d (2,) ^ a (1) _ ^ (0) 

9M TT 87r2 ^ ^ 

M^ S(2-) - -^ E(") (3.4) 

dM 47r2 ^ -^ 



Vertex: 

M— [y(2j) + y(2fe) _^ y(20l ^ ^ y(l) _ i^ y(0) /g gN 

The feature the above equations clearly display, which has already been alluded to in Section 2, 
is the absence of promotion of 1-loop logarithms at two loops for i) the 2-loop vacuum polarization, 
Eq.(3.1), a) the fermion self-energy with a vacuum polarization subdiagram, Eq.(3.4), and Hi) the 
vertex with a vacuum polarization subdiagram, Eq.(3.6). All the rest follow the expected pattern 
of log promotions. Given that the Ward identities provide a relation between Ms and My, but 
not between Mn and anything else, we might then understand the above as a manifest, converse 
statement to the Ward identity, namely, that the renormalization of the vacuum polarization runs 
entirely independently of the renormalization of the self-energy and vertex. Thus, for the 2-loop 
vacuum polarization, although both Iljjif and Iljj^ contain promoted logs of Ms and My (cf. 
Eqs.(2.3) and (2.6)) coming from their respective subdivergences, in the entire 2-loop amplitude 
these promotions cancel. This is as it should be, since otherwise RG equations would imply a relation 
between the mass scales M^ (or My) and Mn, which cannot happen if the renormalization of 
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these amplitudes is to be independent. For the self-energy and vertex with vacuuum polarization 
subdivergences, the same happens: no \n^ M^ occur. The mass derivative of the vertex does give 
a 1-loop vertex, but it is in Landau gauge, and thus finite. For these amplitudes then, we can state 
that there are no genuine 2-loop divergences, and that is a direct consequence of the lack of a Ward 
identity relating the vacuum polarization to any other vertex functions. 

We see then that the renormalization of these amplitudes has, so to speak, split into different 
sectors. A careful study of this structured renormalization pattern should furthermore allow us to make 
predictions about the renormalization at higher loops. In [ll](cf. p. 423), for instance, the statement 
is made that all higher-loop vacuum polarization amplitudes with a single fermion loop (i.e., with 
no lower-loop vacuum polarization subdivergences) do not contain genuine higher loop divergences 
(In^'s, or 1/e^'s in dimensional regularization, etc.). This feature, which may seem fortuitous in other 
renormalization methods, appears naturally in differential renormalization. 

4. Chiral Anomaly 

In this section, we review the computation of the anomalous triangle amplitude {j^{x)ii,{y)j\{z)) 
at one loop [1], and carry it on partially at two loops. Some of the main attractive features of 
differential renormalization are made manifest. We also compare it to another coordinate space, 
regulator- free computation of the anomaly due to K. Johnson [20] for the sake of completeness. 

Since differential renormalization is strictly 4-dimensional and does not introduce any unphysical 
regulator fields, one is able to avoid the complications introduced by standard methods such as 
dimensional regularization (in necessitating ad hoc 75-prescriptions away from d = 4), and Pauli- 
Villars regularization. Furthermore, differential renormalization does not present the anomaly as 
coming from a symmetry that is broken by regularization artifacts, as is usually the case. Instead, we 
shall see that the bare triangle amplitude, without the 7-matrix trace factor in front, has singularities 
which bring in two renormalization scales and these, when the trace is included, lead to a finite 
amplitude with a continuous one-parameter shift freedom given by the quotient of these scales. It 
turns out that this freedom cannot accomodate simultaneously vector an axial conservation. Therefore, 
the triangle amplitude is overconstrained by the two Ward identities. In our scheme, both vector and 
axial symmetries appear manifestly on the same footing, and the shift in the anomaly from one to the 
other is transparent. 
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We begin with the 1-loop computation. The basic elements have already been spelled out in [1]. 
Here we simply sketch the key points. At one loop, the anomalous triangle is (Fig. 3): 



Tf,,,x{x,y) = {j,,{x)ju{y)jUO)) = 2 Tr [757A7a7i.7b7M7c] da— df 



1 



■dr. 



'2/2 {x~yy 
where the factor of 2 reflects the inclusion of the Bose symmetrized diagram. The term 



tabc{x,v) = da^db- 



■a± 



(4.1) 



' y^ {x — y)^ x^ 

is power counting L^^ and thus linearly divergent. This is renormalized in the same way we have 
treated triangles previously, with the difference that now iwo derivatives need to be taken out. This 
has already been done in [1], and we present the final result: 



tabc{x, y) = Fabc{x, y) + Sabc{x, y) ■ 



Fabc is the finite part 
92 



Fabc{x,y) 



dxadyb 
d 
dyt 



a;2y2 (x — yy 



d 

dXa 



hi {^^^^ T°) 0^ 



x^y 



1 



a;2j/2 \^ 4 y (x — j/)2 



1 



r),2y2 



dadbdc - - {S(abdc))a 



1 



(x-yy 

(4.2) 



where () means unnormalized symmetrization in all indices, and Sabc is the renormalization piece: 
1 



Sabc{x,y) = ^TT^ 



^hc r^ Oac ^ 

oxa oyt 



S{x-y)D^-^ 



^ , d d \ / d d \ / a d 

oxa dyaj \oxb dybj \oxc dye 



S{x - y)D 



lnM|a;2 



(4.3) 



We note that, like for the 1-loop fermion self-energy, two different mass scales have been used to 
renormalize the different divergent trace pieces. Not only are we entitled to do that, but in fact one 
more time it will prove crucial that we do so. 

At this point, we apply the traces outside to get Tf^i,\(x,y): 



Tf,i,x{x,y) = R^,^\{x,y) + af,i,x{x,y) 



(4.4) 



where RuiyX comes from the finite part Fabc, and 



4, Ml / d d 



a^ux{x, y) = -167r In — e^^xp ( ^ - ^ ) 5{x)5{y) 



(4.5) 
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It is fairly simple to see that for x ^ y ^ Q, T^i^x is conserved on all channels. However, just 
as in the vertex Wl studied previously, the subtleties are of course in the contact terms. It is also 
w/orthw/hile noting that while Sabc (and thus tabc) indeed contains divergences, the precise 7-structure 
in front arranges things such as to give the final combination In ^ as the only renormalization 
mass dependence, implying that a^uX (and thus T^^a) is actually finite (because M-^ on that 
vanishes). This leads us then to a very clear physical picture: we have a finite anomalous triangle, 
which however contains an ambiguity - in the choice of Ms - coming from a power counting linearly 
divergent Feynman diagram. 

By Fourier transforming into momentum space, one can verify the conservation laws satisfied by 
RfivX on all three channels. Given the form of a^^A found above, the vector and axial Wis on T^i^a 
then read: 

^-T^,x{x,y) = 8tt\1 + 21nMi/M2)e,AMp^— — (5(x)J(y) 

—r^,A(2;,y) = 8^4(1 + 21nMi/M2)e^A.p — ^'5(x)%) (4.6) 

- ig^ + —J T^ux{x,y) = 16^4(1 - 2lnM,/M2)e^,.xpg^—Six)Siy) . 
This is the final and, as it were, a most "democratic" expression of the anomaly: we can tune M1/M2 
so that Tp,yA is conserved either on the vector channels or on the axial one, but never on both; we 
have been able to complete the calculation in a "scheme-free" fashion all the way to the end, without 
having to commit at any point to conservation on a particular channel. 

We now proceed to the 2-loop calculations. The relevant diagrams are indicated in Fig. 4. 
Writing the amplitude as 

T^i^Jxi^^v) - A^,x{x,y)+B^,x{x,y) , (4.7) 

where A^^x contains the contributions of the diagrams with fermion self-energy insertions, and i?p,yA 
contains the contributions of the diagrams with the vertex insertions, we can immediately write, as 
the result of a trivial computation: 



Af,^x{x, y) = ^^ Tr [^5lxlalvlbl,ilc] \da—db , _ n2 ^' 

Oa 5 Ob 7 TT Oc 7: \ Oa :t Ob 7 T7\ Oc ^ 

y^ (y — x)^ x'^ y^ [y — xy x'^ 

To compute the diagrams related to vertex corrections we operate as follows. Let us take, for instance, 

the correction to the vertex at y, which leads to the following integral: 

d'^ud'^v- T^da^^db-. -T^dc-, T^dd-, rrf^e^- (4.9) 

{u — vy w^ [u — yy [y — vy (v — x)^ x'^ 
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Treating the vertex in the standard way, we integrate db by parts over y, and separate a trace and a 
traceless piece, and a surface term. The last two terms do not produce logarithms. It is only the trace 
piece, which is very easy to compute, that brings in a log of the same kind as in Eq.(4.8). Putting 
together the log pieces of the three vertex correction diagrams we get 

B^.x{x,y) = -32^Tr[757A7a7.7fc7M7c] [da^dbj^j—^dc -^ + 

InyHll 1 ^ 1^1;, Hy~xlMl 1 \ 

y^ {y — xy x'^ y^ [y — xy x'^ J 

(4.10) 

where the dots indicate the finite pieces not written explicitly here. Clearly, A^,j\ and B^^a cancel 

exactly except for the fact that the mass scale involved in A^i,\ is A/s whereas the one in B^i,\ is 

My. The one-loop Ward identity shows that the mistmatch is precisely proportional to the 1-loop 

amplitude, thus showing that the entire 2-loop amplitude is finite. The finite parts of B^,^\ are 

extremely lengthy to calculate and we have not found a compact way to perform the computation. 

In fact, the total amplitude is expected to vanish identically since, otherwise, as we now explain, the 

anomaly would get a finite renormalization. 

QED is conformally invariant up to the photon propagator and the appearance of renormalization 
scales. Thus, the anomalous triangle is conformally covariant at one loop because it is finite and 
contains no photon lines. In any coordinate space treatment, this can be verified explicitly to one 
loop, whereas in momentum space calculations, this is obscured because plane waves transform easily 
under translations but not under conformal tranformations. This conformal property, when conjectured 
to all loops, is powerful enough to show the vanishing of all higher-loop triangle amplitudes, in an 
argument due to Baker and Johnson [16]. From the study of conformal transformations on functions 
of three variables, it turns out there is a unique nonlocal, conformal covariant, parity-odd, dimension 
3, VVA tensor [19], and the 1-loop triangle has precisely this form. Now, from uniqueness, if the 
triangle amplitude is covariant to all orders, it means higher-loop contributions also have to have the 
form of the 1-loop triangle. But that would mean a renormalization of the very structure which gives 
rise to the chiral anomaly, and thus a renormalization of the anomaly itself, which is forbidden by 
the Adier-Bardeen theorem. Therefore, all higher-loop contributions to the basic triangle must vanish 
identically. 

Let us finally mention a very nice coordinate space computation of the anomaly which makes no 
use of an ultraviolet regulator, due to Johnson [20]. His starting point is the conformally covariant 

22 



and manifestly finite form of the triangle gotten by acting out the 7-matrix trace in Eq.(4.1) from 
the beginning. Because there is, apart from this nonlocal structure, a unique contact term with 
precisely the same dimension and conformal properties [19], the triangle amplitude can have this local 
ambiguity, and its coefficient can be chosen to give conservation on either channel but not both. If 
we want conservation on the x channel, for instance, we have: 



d 



T,.x(x, y) + ae^"^^- j^-T- _ __ ) S{x)S{y) 



d_ d_ 

^ dx'^ dy" 



0, (4.11) 



where we have shown explicitly the contact term in question, with coefficient a to be determined. 
To find a, one then integrates d^T^i^x in x and y against the "test" function x"y^ . This integral, 
remarkably, is entirely determined from the long-distance behavior of T^^a, and the only cutoff needed 
to perform it is an infrared one. Johnson also points out that this situation is reminiscent of the 
Poisson equation in classical electrostatics, where the coefficient of a contact term in a differential 
equation (the charge of a point particle) is also determined by the long-distance behavior of a field 
(Gauss' theorem). 

5. Conclusion 

In this paper, we have presented a detailed analysis of the differential renormalization of QED up 
to two loops. All computations are reasonably simple when compared to other approaches. 

Our results can be summarized as follows. We have found explicit expressions for the 1- and 2- 
loop renormalized amplitudes of the vacuum polarization, self-energy of the fermion and vertex. Only 
some finite parts of the latter are left in the form of integrals. Ward Identities are verified and provide 
independent checks of the computations. Furthermore, the amplitudes obey renormalization group 
equations which yield the beta function and the various anomalous dimensions of all the basic fields in 
the theory. This renormalization group equations display a natural organization due to the freedom to 
use different renormalization subtractions for quantities not related by Wis. We have also presented 
the treatment of the anomalous triangle which, in the very spirit of differential renormalization, is 
regulator-free. It sides with other techniques which are based on the interplay between potential 
counterterms and symmetries rather than on explicit computations on a regulated theory. 

One is rapidly enticed by the ease of computation of differential renormalization. It is our opinion 
that the method also deserves further consideration because of its natural treatment of chiral problems. 
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Appendix A. Basic Differential Renormalization Identities 



The following are the DR identities used to renormalize all of the amplitudes presented in the 
text: 



-r = -Tn- 



1 Inx'^M'^ 



1 



1 Inx^M^ 
rDD ^ 



x^ 32 X- 

= — n — ^^ 

a;4 8 a;2 



(A.l) 
(A.2) 
(A.3) 



Appendix B. Two-loop renormalized vertices 



We present here the final expressions for the 2-loop renormalized vertices. In the expressions 
below, a derivative w.r.t. z^ means -^^ = — -jAr — ttti- 



'a\ 5/2 



IGtt'^ Vtt/ 



(7b7M7a) 



Inw^M^ ^ 1 
-rl - 



y\x -yY \ 4 



dz^ y-^{x — yY x^ 

•2 In^y^^l + 21ny2M2 



^5ab5'^'^\x)n 



y" 



(B.i) 



T/p'^)(x,y) 



967r3 
d 



/a\5/2 f 



dy'' \x'^y 



1 



dx'^ 



da{dpd„ - 5p„n)L{x - y) 



9/1 1 

^db^idpdo- - SpaD)L{x -y)] + 



16 



la dadf. 



ln(x - yyMfi 

{x — J/)2 



(B.2) 
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where L{x) = Inx^^^ _ i ^^^2 2,2^^/2 



V!l''\x,y) 



3^ ©'^V'^'"^^ {(7.7.7.7c7a) ^ 



1 



.5f 



1 1 



rda- 



1 



r^d- 



1 



, (x - yf ln(a; - y)'^M'^ 



1 



(a: — y)'^v'^ v? {u — v)^ {x — u)^ (v — y)^ 



(X — J/j'^ w*^(m — Wj^ 



16^ (- j (7<.7.7a) I ^ 



r^a 



(a; — yY y'^x'^ 
647r Vtt 



0''%,^(^)(x): 



(a;-y) 



1 ,^ ^ &ad. 1 Iny^M^ 
^{dadd - -rn)-^ , 



x" y" 



(B.3) 



— i /a\5/2 



647r'^ Vtt 

- 9a ^r 



(7p7a7A*767d7p7c) j dudv \-^ 



da^^db 



1 



r^e 



1 



1 



1 



db — {dcdd 1-0)7 — +dd- 



x^ u^ [x — w)2 (m — w)2 (u — y)2 [u — yY 
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(u — yY {u — y) 



1 ,r, r. ^Cb„. 1 



(a; — vY{u — vY{v — yY 
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167r^ Vtt 
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11 ln(a; - yYM"^ 
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x^ y"^ (x — yY 
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In (x - yyM{, + 2 ln(x - yYM 
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Figure captions: 

Figure 1: 1-loop IPI diagrams of QED. 

Figure 2: 2-loop diagrams of QED. 

Figure 3: 1-loop anomalous triangle diagram. 

Figure 4: 2-loop contributions to anomalous ABJ amplitude. 
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